Universal characteristics of resonant-tunneling field emission from nanostructured surfaces 
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We have performed theoretical and experimental studies of field emission from nanostructured semiconductor 
cathodes. Resonant tunneling through electric-field-induced interface bound states is found to strongly affect 
the field-emission characteristics. Our analytical theory predicts power-law and Lorentzian-shaped current- 
voltage curves for resonant-tunneling field emission from three-dimensional substrates and two-dimensional 
accumulation layers, respectively. These predicted line shapes are observed in field emission characteristics 
from self-assembled silicon nanostructures. A simple model describes formation of an accumulation layer and 
of the resonant level in these systems. 
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I. INTRODUCTION 

Field emission describes the ejection of electrons from 
a conducting substrate into vacuum via electron tunneling 
through the surface potential barrier.'-^ As tunneling is a sen- 
sitive probe of the electronic spectrum, field emission has 
often been applied as a useful spectroscopic tool for bulk- 
material characterization.^^ Experiments using nanostructured 
emitters'*'^ reveal deviations from conventional theory sig- 
nifying that, in these systems, tunneling occurs through dis- 
crete energy levels formed in the nanotips due to spatial con- 
finement. In addition to revealing intriguing properties of 
electron transport at the nanoscale, the intense and highly 
coherent electron beams extracted from such nanotips have 
found useful applications in projection electron holography.^ 
More recently, a number of field-emission studies of 2D 
(quantum film), ID (quantum wire),^'^ and OD (quantum dot)'' 
semiconductor systems have also demonstrated clear devi- 
ations from Fowler-Nordheim behavior In particular, dis- 
crete current peaks, an anomalously low threshold field for 
electron emission and regions of negative differential con- 
ductance have all been observed in the current-voltage (/-V) 
field-emission characteristics. All of these features are as- 
sociated with the formation of quantized bound states in the 
low-dimensional semiconductor cathodes. These experimen- 
tal results have stimulated a number of theoretical studies of 
electron field emission from and through quantum-confined 
electronic states. '"'' '''^-'^''^ 

In this work, we consider a generic system consisting of a 
bulk semiconductor substrate (S) with a nanostructured sur- 
face where nanometer-sized tips ('whiskers') have formed. 
In addition to the energy baiTier between electronic states in 
the nanostructure (N) and vacuum (V), quantum confinement 
gives rise to a further barrier at the substrate-nanostructure 
(S-N) interface. Electric-field penetration into the semicon- 
ductor material leads to band bending'^ '^ and the formation 
of 2D bound states at both interfaces. The evolution of bound- 
state formation is illustrated in Fig. [T] for a set of model pa- 
rameters. We focus on the situation shown where, initially, a 
2D bound state is formed at the S-N interface. See Fig. [Ha). 
As its energy lies below the substrate conduction-band edge, 
no resonant tunneling can occur via this state. As the electric 
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FIG. 1: Energy band diagram of a realistic model for a vacuum- 
nanostructure-substrate system based on a silicon nanostructure 8 
nm high and 10 nm in diameter. The external electric field F is 
increasing in successive panels (a)-(c). V, N, and S correspond to 
the vacuum, nanostructure and substrate regions, while and iJnv 
correspond to the lowest energy bound states located at the substrate- 
nanostructure and nanostructure-vacuum interfaces respectively. 



field is increased further, quantization of electron energy also 
occurs within the potential well formed at the N-V interface, 
as shown in Fig. [Hb). The corresponding bound states can 
mediate resonant tunnelling both from the 3D substrate elec- 
tron system and also from the 2D bound states formed in the 
S-N quantum well. Here we present results from our analyti- 
cal theory which describes resonant-tunneling field emission 
from such 3D and 2D systems. These studies provide new in- 
sight into universal features exhibited in field emission from 
nanostructured surfaces and which are confirmed by our mea- 
surements of field-emission currents from self-assembled sil- 
icon nanostructures. In particular, resonant tunneling from a 
3D substrate via a 2D interface bound state is found to give 
rise to a power law / oc V^^^ in the ascending part of the res- 
onance peak in the I-V curve. The clear contrast to Fowler- 
Nordheim theory, which predicts / oc y^, allows for unam- 
biguous identification of the 3D-vacuum resonant-tunneling 
mechanism. Alternatively, when resonant tunneling through 
the 2D interface bound state occurs from a 2D accumulation 
layer, the peak exhibits a Lorentzian line shape. We start by 
describing the theoretical method and results. This is followed 



2 



by details of our sample preparation, measurement setup, and 
the obtained field-emission data. 



11. THEORETICAL ANALYSIS 

The general expression for the field-emission current den- 
sity from a cathode as a function of temperature T and voltage 
V is given by 
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Here T{E, V) denotes the energy-dependent quantum- 
mechanical transmission function from electrode to vacuum, 
and N{E, T, V) is the supply function of the electron reser- 
voir, i.e., the density of electrons in transverse plane-wave 
states having energy E in their motion perpendicular to the in- 
terfaces. The substrate-vacuum voltage drop V is proportional 
to the applied electric field F, with the constant of proportion- 
ality being dependent on the specific experimental geometry 
and system under consideration. We proceed now to discuss 
the form of the transmission and supply functions employed 
in our transport calculations before stating the results of our 
analytical theory. 

Resonant tunneling through a 2D bound state, e.g., the one 
formed at the N-V interface shown in Fig.lT] can be described 
phenomenologically'^ by a Lorentzian energy dependence of 
the transmission coefficient. 
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where Env and Fnv are the resonance energy and life-time 
broadening, respectively, of the bound state at the N-V in- 
terface. Tpk is the transmission probability at resonance. As 
the double-barrier system under consideration is highly asym- 
metric, with a high (low) barrier at the N-V (S-N) interface, 
the peak transmission is dominated by Fowler-Nordheim tun- 
nelingi^ through the N-V barrier The relation 



(-J -^pk Fnv ~ -^nv -^nv(-^nv) 



can then be derived,'^ with the transmission function for tun- 
neling through the N-V barrier denoted by T„y{E). 

The explicit form of the supply function depends on the 
dimensionality of the electron reservoir For field emission 
from a 3D substrate, it is given by the familiar expressior^ii 
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with fee denoting the Boltzmann constant, m being the ef- 
fective mass of electrons in the substrate, and /is the equilib- 
rium chemical potential of the substrate. The supply function 
is modified, however, when a 2D system such as the accu- 
mulation layer formed at the S-N interface acts as the elec- 
tron reservoir While the energy of electron motion trans- 
verse to interfaces is quantized in a quantum well, coupling 
of the bound state to substrate and nanostructure lead to a fi- 
nite life time of electrons in the accumulation layer. Hence it 
supplies electrons with energies in a narrow range Fsn around 
the bound-state energy E^n with a Lorentzian distribution,— 
yielding 
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Having specified transmission and supply functions, we are 
in the position to calculate the field-emission characteristics 
using Eq. ([T]i. To obtain our analytical results, we make the 
following assumptions: (i) presence of a highly degenerate 
electron gas in the reservoir, which means that we use the 
T ^ limiting expression for the supply functions, and 
(ii) triangular shape of the N-V quantum well, which allows 
us to express E„y in terms of the applied electric field F as'i 
Eny = {h'^TT^/2me^y/^{eF)^/^. We find then the following 
expressions for the I-V characteristics of resonant-tunneling 
field emission from a 3D substrate and from a 2D accumula- 
tion layer, respectively; 



.?rtfe(^) 



(eV) 3 exp 



Vo 
V 



' (Ms - E„,) 9 {E,, - Ei,,) 



(6a) 



r 



iRTFE(^) ~ iRTFE(^) 



(3D) 



Fsn "l~ Fn 



(Ais - E,,) 



1 



(6b) 



Here k denotes the dielectric constant of the substrate ma- 
terial, D is the distance over which the cathode potential 
is dropped, and Vb is a voltage scale related to the effec- 



tive work function of the vacuum barrier, which is modi- 
fied due to energy quantisation. The Heaviside step func- 
tions 9{x) in Eqs. (|6]l enforce proper phase-space conditions 
for resonant-tunneling transport to occur, e.g., cutting the 
current off sharply when the (in general voltage-dependent) 
resonant-level energy E^^ falls below the conduction-band 
bottom ii'bs- The conspicuous difference in the voltage- 
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FIG. 2: 1.5 /imxl.5/im AFM surface plot of the Si nanostmctures 
formed by electron beam annealing at 1000 °C for 15 s. The rate at 
which the sample was heated and cooled was 5 °C/s. 

dependent pre-exponential factor of Eq. (l6a] l compared to the 
Fowler-Nordheim formula [{eV)^^^ as opposed to (eV^)^] 
arises from the field dependence of the bound-state energy 
Eav oc (V/D)"^^^ formed in the triangular N-V quantum well. 

III. EXPERIMENTAL PROCEDURE 

We now proceed to describe details of our experimen- 
tal study which examines the field-emission properties of a 
self-assembled silicon nanostructure formed on a silicon sub- 
strate.''^ Detailed analysis of the I-V field-emission charac- 
teristics indicates the presence of two classes of current peaks 
which occur at different regions of electric field. Based on 
our understanding gained from theoretical analysis described 
above, we can relate these two classes of current peaks to res- 
onant tunnelling from the 3D substrate and from 2D bound 
states located at the N-S interface, respectively. In both cases, 
tunnelling is mediated by a second set of 2D bound states lo- 
cated at the N-V interface. 

Substrates used in the present study were n-type Si (100) 
wafers (P-doped, conductance 1-10 flcm , corresponding to 
a dopant density of lO^^-lO'^^ cm^'^). Nanostructure growth 
results from annealing of the untreated substrates to 1000 °C 
for 15 s using a 20 keV raster-scanned electron beam. Details 
of the fabrication procedure and mechanism for nanostructure 
growth have been reported elsewhere.^" A characteristic AFM 
image of the nanostructured Si substrate surface is shown in 
Fig. 121 The self-assembled nanostructures are square-based 
and distributed randomly over the surface with a density of 1 1 
/im~^. High-resolution AFM and TEM studies^ indicate an 
average nanostructure height of 8 nm and a base length distri- 
bution in the range 8-60 nm. The relevance of these structures 
to the study of 3D and 2D resonant-tunnelling field emis- 
sion phenomena of the type discussed theoretically above, be- 
comes apparent when one considers an energy band diagram 
of the S-N-V system. In the first case, the S-N interface was 
modelled by an offset in the conduction band equal to the con- 
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FIG. 3: I-V field emission characteristics for Si nanostructures 
formed on n-type Si (100) substrate. The three curves correspond to 
successive conditioning cycles, where 1 corresponds to the first and 
3 the final conditioning cycle . The solid lines that overlay curves 
2 and 3 correspond to a fit to Fowler-Nordheim theory. For clar- 
ity, each curve has been offset by 1 jiA . The inset shows Fowler- 
Nordheim plots of curves 1 and 2. The arrows indicate the positions 
of the five most prominent peaks, while the solid lines indicate the 
best fit for the experimental data in the Fowler-Nordheim emission 
region. Curve 2 has been offset by 1 along the y-axis. 



finement energy of the lowest-lying transverse bound state 
in the nanostructure. Self-consistent solutions of the Pois- 
son and Schrodinger equations were then used to calculate the 
conduction-band-bottom profile and to find regions of elec- 
tron energy quantization.-^ Given the size distribution of the 
nanostructures, and assuming a rectangular quantum well, the 
offset in the conduction band was calculated to be in the range 
0.6-31 meV with bound states located at the S-N interface 
forming only in those structures smaller than 12 nm. Fig.[T] 
shows the results of our simulations for a 10 nm diameter sil- 
icon nanostructure. 



IV. RESULTS AND DISCUSSION 

Characterisation of the field-emission properties of these 
self-assembled nanostructures was performed at room temper- 
ature in a parallel-plate diode configuration. Electrical isola- 
tion and cathode-anode separation was achieved using a 100 
/im thick PTFE film, and measurements were performed un- 
der vacuum of 1 x 10~^mbar It was found necessary to run 
several emission cycles at high emission cuiTents (5 /iA ) in 
order for the I-V characteristics to become stable and repro- 
ducible. The measured anode voltage can be related to the 
voltage drop V, defined in the equations above, via a lever- 
arm relationship. The I-V curves plotted in Fig. [3] show the 
change in emission current during the course of a condition- 
ing cycle. The nanostructured substrate was maintained un- 
der vacuum between successive cycles. Following incomplete 
conditioning (curve 1) the emission characteristics are typi- 
fied by two regions of exponentially rising current (labelled 
a and c, respectively), separated by a plateau in emission cur- 
rent (region b). With continued conditioning however, distinct 
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FIG. 4: (a) Field emission resonance region following subtraction 
of the background Fowler-Nordheim emission current as fitted to the 
I-V characteristics displayed in (b). (c) Semi-logarithmic plot of 
experimental data shown in (b). Solid lines indicate the best fits to 
current peaks FN^ and FN^. 

current peaks develop within the current plateau region (see 
curves 2 and 3 of Fig. |3]l. The arrows indicate the positions 
of the five most prominent peaks. While the peak positions 
and relative intensities were found to be repeatable for each 
sample, a significant deviation in the number, position and in- 
tensity of peaks was observed between samples. While it is 
not possible to identify the individual nanostructure involved 
in electron emission, AFM images of the substrate surface, 
taken before and after field emission measurements, appear 
identical. Further, following exposure to air, it was found nec- 
essary to again condition the cathodes in order for the current 
peaks to develop, suggesting that conditioning is related to re- 
moval of adsorbates from the cathode surface. 

The current peaks are also clearly observed when plotting 
the I-V characteristics on a Fowler-Nordheim plot [ln(//y^) 
versus 1/V~\, as shown in the inset of Fig. |3] for curves 1 
and 2. The Fowler-Nordheim plot also reveals highly lin- 
ear behaviour at anode voltages > 550 V, suggesting electron 
emission through conventional (direct, non-resonant) Fowler- 
Nordheim tunnelling. The solid lines in Fig. |3] which dis- 
play the Fowler-Nordheim relation fitted to the experimental 
data in this high voltage region, demonstrate excellent agree- 
ment between the model and observed emission current. Cur- 
rent peaks in the I-V characteristics of substrates containing 
multiple cathodes could be related to the activation and de- 
activation of additional emission sites. However, the leading 
edge of each peak marked in the inset of Fig. [3] is found to 
be highly non-linear, suggesting conduction through a mech- 
anism other than conventional Fowler-Nordheim tunnelling, 
in contrast to what would be expected from the activation of 
additional emission sites. This is in agreement with measure- 



ments performed using a phosphor-coated anode, which sug- 
gests electron emission from a single emission site over the 
range of anode voltage considered. 

In order to study the current peaks in more detail, the back- 
ground current fitted to the high-voltage Fowler-Nordheim 
tunnelling region was subtracted from the experimental data, 
as shown in Fig. |4|a). The corresponding I-V characteris- 
tics are displayed in Fig. |4|b). Four current peaks are ob- 
served (labelled 1^) which supply an emission current 1.8 
times greater than that predicted from conventional Fowler- 
Nordheim tunnelling. In addition to these symmetric current 
peaks, we find two highly asymmetric peaks, labelled FN^ 
and FN^ in Fig. HI The presence of the asymmetric peaks be- 
comes particularly apparent when the I-V characteristics of 
Fig. IHb) are displayed on a semi-logarithmic plot, as shown 
in Fig.|3c), where discontinuities aligned with peaks FN^ and 
FN'^ are clearly observed. It is possible to match the experi- 
mentally observed field-emission resonance peaks to our the- 
oretical results presented earlier Fits to the line shape of peaks 
FN^ and FN^, shown by the solid lines in Fig.|4|c), yield pre- 
exponential factors of l/^'^** and 1/1-704^ respectively. This 
is very close to the theoretically expected power law V^/^ 
predicted for resonant tunneling from the 3D substrate states. 
Furthermore, the Lorentzian line shape of resonance peaks la- 
belled 1^ [see Fig.Ua)] is in good agreement to our theoreti- 
cally predicted lineshape corresponding to resonant tunneling 
from the S-N accumulation layer The excellent agreement 
between experiment and theory is a strong indication that 
resonant-tunneling transport as illustrated in Fig. [T] indeed 
occurs in our samples^"*. Additional numerical calculations 
based on previously described methods' "'^'''^'^^'^"^ and further 
modelling of sample electrostatics could yield the bound-state 
energies and their life-time broadening but these extended 
studies are beyond the focus of this work. 

In conclusion, we have studied field emission from nanos- 
tructured surfaces that is strongly affected by resonant tun- 
neling from both a 3D substrate and from 2D accumulation 
layers. Universal features in the current-voltage charateris- 
tics are derived theoretically and confirmed experimentally, 
and allow for identification of resonant-tunneling peaks with 
specific microscopic tunneling mechanisms. While the exper- 
imental part of this study considers resonant tunnelling from 
silicon nanostructure field-emission cathodes, the analytical 
expressions derived are quite generic and are applicable to a 
range of semiconductor structures used in the study of field- 
emission phenomena. 

This work was performed under research contracts to the 
New Zealand Foundation for Research, Science and Tech- 
nology (C05X0008, 'Advanced industrial materials at the 
nanoscale'). 
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